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Abstract 

The dynamics of symmetry breaking during out of equilibrium phase tran- 
sitions is a topic of great importance in many disciplines, from condensed 
matter to particle physics and early Universe cosmology with definite experi- 
mental impact. In these notes we provide a summary of the relevant aspects 
of the dynamics of symmetry breaking in many different fields with empha- 
sis on the experimental realizations. In condensed matter we address the 
dynamics of phase ordering, the emergence of condensates, coarsening and 
dynamical scaling. In QCD the possibility of disoriented chiral condensates 
of pions emerging during a strongly out of equilibrium phase transition is dis- 
cussed. We elaborate on the dynamics of phase ordering in phase transitions 
in the Early Universe, in particular the emergence of condensates and scaling 
in FRW cosmologies. We mention some experimental efforts in different fields 
that study this wide ranging phenomena and offer a quantitative theoretical 
description both at the phenomenological level in condensed matter, intro- 
ducing the scaling hypothesis as well as at a microscopic level in quantum 
field theories. The emergence of semiclassical condensates and a dynamical 
length scale is shown in detail, in quantum field theory this length scale is 
constrained by causality. 

The large N limit provides a natural bridge to compare the solutions in 
different settings and to establish similarities and differences. 
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I. PHASE ORDERING DYNAMICS: AN INTERDISCIPLINARY FASCINATING 

PROBLEM 



The dynamics of non-equilibrium phase transitions and the ordering process that occurs 
until the system reaches a broken symmetry equilibrium state play an important role in 
many different areas. In condensed matter physics binary fluids, ferromagnets, superfluids, 
and liquid crystals, to name a few, are examples of systems in which the dynamics of phase 
transitions out of equilibrium are studied experimentally. 

Experiments in these systems have provided a solid basis for the description of the 
dynamics of phase ordering: in binary fluids or alloys upon a sudden temperature drop 
below the critical temperature, the two fluids begin to separate, regions of different fluid 
concentrations are separated by domain walls. In superfluids, rapid cooling leads to a 
network of vortices and in liquid crystals to many different topological defects. 

In cosmology defects are conjectured to be produced during Grand Unified Theory (GUT) 
or the Electro-weak (EW) phase transition can act as seeds for the formation of large scale 
structure and the dynamics of phase ordering and formation of ordered regions is at the heart 
of Kibble's mechanism of defect formation |1]-§J. Current and future measurements of Cosmic 
Microwave Background anisotropies could provide distinct evidence for these phase transi- 
tions whose dynamics could have influenced structure formation Q]. At even lower energies, 
available with current and forthcoming accelerators, the Relativistic Heavy Ion Collider 
(RHIC) at Brookhaven and the Large Hadron Collider (LHC) at Cern the phase transitions 
predicted by the theory of strong interactions, Quantum Chromodynamics (QCD) could 
occur out of equilibrium resulting in the formation of coherent condensates of low energy 
Pions. These conjectured configurations known as 'Disoriented Chiral Condensates' are sim- 
ilar to the defects expected in liquid crystals or ferromagnets in condensed matter systems 
and their charge distribution could be an experimental telltale of the chiral phase transition 
of QCD ||. Whereas the GUT phase transition took place when the Universe was about 
1CT 35 seconds old and the temperature about lO 23 ^, and the EW phase transition occured 
when the Universe was 10~ 12 seconds old and with a temperature 10 15 i^, the QCD phase 
transition took place at about 10 -5 seconds after the Big Bang, when the temperature was a 
mere lO 12 ^. This temperature range will be probed at RHIC and LHC within the next very 
few years. The basic problem of describing the process of phase ordering, the competition 
between different broken symmetry states and the formation and evolution of condensates 
on the way towards reaching equilibrium is common to all of these situations and fields. The 
tools, however, are necessarily very different: whereas ferromagnets, binary fluids or alloys 
etc, can be described via a phenomenological (stochastic) classical description, certainly in 
quantum field theory a microscopic formulation must be provided. In these lectures we 
describe a program to include ideas from condensed matter to the realm of quantum field 
theory, to describe the non-equilibrium dynamics of symmetry breaking and the process of 
phase separation and phase ordering on a range of time and spatial scales of unprecedented 
resolution (for example in QCD the time scales < 10~ 23 seconds, spatial scales < 10~ 15 
meters, in cosmology the time scales are of order 10~ 32 seconds and spatial scales smaller 
than 10~ 44 meters) that require a full quantum field theoretical description. 

We begin the excursion into these timely fields by first providing a brief quantitative de- 
scription of the relevant setting and whenever possible the experimental situation associated 
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with them in three main areas: Condensed Matter, Ultrarelativistic Heavy Ion Collisions 
and Early Universe Cosmology. This quantitative discussion will be followed by a more 
qualitative description of some of the main theoretical ideas, techniques and tools. 

A. Condensed Matter: 

A description of phase transitions and critical phenomena in equilibrium begins by rec- 
ognizing an order parameter which is a thermodynamic ensemble average of a macroscopic 
variable that determines the different macroscopic states of the system. For example in fer- 
romagnets the order parameter is the average magnetization, above a critical temperature 
it vanishes and it is non-zero below the critical temperature, in superfluids is the conden- 
sate density, in superconductors the density of Cooper pairs, etc. |J. Phase transitions in 
equilibrium are fairly well understood and described by the theory of critical phenomena 
U which combined with the renormalization group provides a very succesful description of 
phase transitions. The theory of critical phenomena and the renormalization group pro- 
vide a very robust description of universality classes: many systems that are very different 
behave similarly near critical points, these universality classes are divided by for example 
the dimensionality of the order parameter, the dimensionality of space, and the symmetries 
of the underlying microscopic Hamiltonian. An important concept in critical phenomena 
is the correlation length, take for example a spin one-half ferromagnet, the microscopic 
Hamiltonian has an up-down symmetry, the energy remains the same if all spins are flipped. 
Focus at a particular point of the sample where the spin is up. The correlation length is 
the distance over which the spins are correlated, i.e. the distance from this up-spin over 
which the neighboring spins are also up. As the critical temperature is reached from above 
this correlation length grows reaching a macroscopic size (diverging) at the critical temper- 
ature. As the system is cooled below the critical temperature a phase transition occurs: 
there appears a net overall magnetization and at low temperatures all spins are either up 
or down, the up-down symmetry is spontaneously broken |J. This phase transition occurs 
in equilibrium when the microscopic relaxation time scales are shorter than the time scale 
of cooling the system, thus at all times the system is in local thermodynamic equilibrium. 
At very high temperatures typically the disordered phase prevails, all spins are oriented 
at random and the average magnetization vanishes. As the critical temperature is reached 
regions of correlated spins appear and become of macroscopic size as the correlation length 
diverges and the spin system begins to order. In this region the thermodynamic quantities 
become insensitive to the short distance details such as crystalline lattices, lattice spacing 
and the nature and strength of the interaction between spins as the physics is determined 
by the correlation of spins over large distances. Near the critical point, the short distance 
length scales are irrelevant for macroscopic phenomena and long-wavelength physics is com- 
pletely determined by the correlation length £(T). Macroscopic thermodynamic quantities 
and susceptibilities near the critical temperature only depend on the length scale £(T). 

This is the basis of the static scaling hypothesis which is confirmed experimentally in a 
wide variety of systems and is theoretically supported by the renormalization group approach 
to critical phenomena M . The static critical phenomena associated with second order phase 
transitions that occur in local thermodynamic equilibrium is fairly well understood via the 
renormalization group (and other alternative approaches) |J. 
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Consider the alternative scenario in which a ferromagnet is held at very high temperature 
in the disordered phase and suddenly it is cooled below the critical temperature on time 
scales shorter than those associated with relaxational phenomena. Now the spin system 
must evolve towards the ordered phase far away from equilibrium. Unlike the case of static 
(local thermodynamic equilibrium) critical phenomena, the case of out of equilibrium phase 
transitions require a novel set of ideas and tools to describe the dynamics of the process of 
phase ordering. 

There is now a large body of theoretical and experimental work in phase ordering dy- 
namics in condensed matter systems 0- [IU|. Although ultimately the tools to study similar 



questions in quantum mechanical many body systems will be different, the main physical 
features to describe are basically the same: as the system cools down suddenly below the 
critical temperature correlated regions (of spins in a ferromagnet or of condensate in a Bose 
superfluid) begin to form. These correlated regions are separated by 'walls' or other struc- 
tures. Inside these regions an ordered phase exists which eventually grows in time to become 
macroscopic in size. 

Before attempting to describe the manner in which a given system orders after being 
cooled through a phase transition an understanding of the relevant time scales is required. 
Two important time scales determine if the transition occurs in or out of equilibrium: the re- 
laxation time of long wavelength fluctuations (since these are the ones that order) r re i{k) and 
the inverse of the cooling rate t coo i = T(t)/T(t). If r re i(k) « t coo i then these wavelengths 
are in local thermodynamical equilibrium (LTE), but if r re i(k) » t CO oi these wavelengths 
fall out of LTE and freeze out, for these the phase transition occurs in a quenched manner. 
These modes do not have time to adjust locally to the temperature change and for them the 
transition from a high temperature phase to a low temperature one occur instantaneously. 
This description was presented by Zurek [ITJ analysing the emergence of defect networks 
after a quenched phase transition. Whereas the short wavelength modes are rapidly ther- 
malized (typically by collisions) the long- wavelength modes with k << 1/£(T) with £(T) the 
correlation length (in the disordered phase) become critically slowed down i.e. their relax- 
ation time becomes extremely long near the critical point. As T — > T c + the long wavelength 
modes relax very slowly, they fall out of LTE and any finite cooling rate causes them to 
undergo a "quenched" non-equilibrium phase transition. As the system is quenched from 
T > T c (disordered phase) to T « T c (ordered phase) ordering does not occur instanta- 
neously. The length scale of the ordered regions grows in time (after some initial transients) 
as the different broken symmetry phases compete to select the final equilibrium state. A dy- 
namical length scale typically emerges which is interpreted as the size of the correlated 
regions, this dynamical correlation length grows in time to become macroscopically large 
|7HT0|. Just as in static critical phenomena, the emergence of this dynamical correlation 
length leads to the dynamical scaling hypothesis, that the approach to equilibrium and the 
kinetics of phase ordering is solely determined by this length scale. 

Experiments in binary fluids for example, study the growth of these correlated regions 
by light scattering |T2j much in the same manner as the onset of ferromagnetism is studied 
via neutron scattering. The growth of the domains is characterized by the dynamical length 
scale. As a function of time the scattering of light becomes stronger for longer wavelengths 
i.e. smaller wave- vectors, until eventually at very long times a Bragg peak at zero momentum 
emerges signaling the macroscopic ordering of the system. This growth of domain structures 
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during the dynamical process of phase ordering is referred to as "coarsening" |7|-[10|. This 



mechanism with a clear experimental realization in condensed matter is at the heart of the 
Kibble-Zurek p|,p|JTl|j3|jr^] scenario of the dynamics of symmetry breaking in cosmological 
phase transitions [|l]J^|n]|||n| . In this scenario a "network" of defects emerges after a phase 
transition that occurred strongly out of equilibrium with a density of about one defect per 
initial correlation length. Computer simulations reveal that this network evolves and a 
scaling regime emerges |fl. 

This idea had sparked an intense effort to reproduce "cosmological phase transitions 
in the laboratory". The dynamics of phase ordering had been studied in liquid crystals 
whose symmetry group is rather similar to that of particle physics models. The experiments 
produced non-equilibrium phase transitions both by suddenly varying the pressure and the 
temperature (pressure and temperature quenches) p^JT^|] and confirmed at least in a qualita- 
tive manner the main features described by this scenario of dynamics of symmetry breaking. 
More recently a new set of experiments had sought to provide a more detailed picture of the 
dynamics of symmetry breaking phase transitions out of equilibrium and to simulate in the 
laboratory what is thought to be the situation in cosmological phase transitions. Original 
experiments focused on studying the dynamics of phase ordering after a pressure quench 
in superfluid 4 He |]TB[ by measuring second sound (i.e. entropy disturbances) which only 
propagate in the superfluid component (the broken symmetry phase). The interpretation of 
results in these experiments were overshadowed by induced turbulence during the quench 
and spurious phase separation due to imperfections of the walls. More recently a new set of 
experiments were carried out that seem to lead to cleaner interpretations. 

In these ingenious experiments fIT|]l8|| a small sample of superfluid 3 He, whose order 



parameter has a group structure very similar to some particle physics models, was heated 
locally by neutron irradiation via the nuclear reaction 

3 He + n = 3 He~ + p + 764 keV 



the energy released heats a small portion of the liquid Helium into the normal state and rapid 
diffusion of the quasiparticles cools this region back into the superfluid phase very rapidly, 
thus providing a quench from a normal (disordered) phase into the superfluid (ordered) 
phase. The resulting domain structure is then studied via NMR and a qualitative agreement 
with the picture of the symmetry breaking dynamics seem to emerge from these experiments. 
Thus these beautiful experiments in condensed matter reproduce in a qualitative manner 
the physics of a "little Big Bang" and provide controlled experimental framework to test the 
concepts associated with the dynamics of symmetry breaking. 

These ideas of the emergence of correlated regions that grow in time and become macro- 
scopic during non-equilibrium phase transitions has been recently invoked as a potential 
signature of the chiral phase transition in QCD, the theory of the strong interactions. 



B. Chiral symmetry breaking in QCD and "disoriented chiral condensates" 

Quantum Chromodynamics (QCD) is the theory of strong interactions, with the funda- 
mental degrees of freedom being the quarks and gluons. Quarks, however are confined inside 
hadrons and do not exist as individual, isolated particles in vacuum. However there is now 
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a wealth of theoretical evidence including very convincing lattice results that indicate that 
at temperatures above T « 150Mev quarks and gluons become free and form a quark-gluon 
plasma. The lattice results are supported qualitatively and quantitatively by phenomenolog- 
ical models fl9j . In fact the numerical evidence supports the picture of two phase transitions: 
the quark-gluon plasma or confining-deconfining phase transition in which the quarks and 
gluons become confined into hadrons and the chiral phase transition that leads to the low 
energy description in terms of pions. The low energy limit of QCD is dominated by the light- 
est up and down quarks u, d with masses m u ~ 5Mev; « 7 — lOMev. These mass scales 
are much smaller than the natural scale of QCD, Aqcd ~ lOOMev at which QCD becomes 
strongly coupled. In the limit of vanishing up and down quark masses, the QCD Hamilto- 
nian possesses a global chiral symmetry corresponding to rotating independently the right 
and left handed components of the spinors that describe the quark fields. This symmetry is 
SU (2) i g) SU{2) R and in the low energy world is spontaneously broken down to SU (2) L+R 
with the charged and neutral pion isospin triplet being the (quasi) Goldstone bosons associ- 
ated with the breakdown of this symmetry. The small mass of the pions (~ 135Mev), on the 
hadronic scale is a result of the small mass of the up and down quarks on the QCD scale, 
which breaks explicitly chiral symmetry. This is the chiral phase transition. The lattice 
results seem to indicate that the two transitions, deconfining and chiral symmetry breaking 



are very close in temperature and may actually happen at the same temperature [19 



Whereas the deconfining phase transition does not seem to be characterized by a natural 
order parameter, the chiral transition is described by the non-vanishing of the chiral con- 
densate < qq > with q = (it, d) with < • • • > refering to the vacuum expectation value or the 
thermodynamic ensemble average at finite temperature. Although this transition(s) have 
taken place when the Universe was at a temperature of 150Mev about 10~ 5 seconds after 
the Big Bang, the Relativistic Heavy Ion Collider (RHIC) at Brookhaven to start operation 
at the end of 1999 and the Large Hadron Collider at Cern (around 2004) will probe this 
transitions by colliding heavy ions. 

RHIC will accelerate and collide from protons up to 250 Gev and ions of up to the heaviest 
nuclei with collision energies of about 100 Gev per nucleon for Au nuclei. The phenomenon 
of nuclear transparency observed in nucleon-nucleon collisions leads to the conclusion that 
about half the energy of the collision is carried away by the nuclei and about half the energy 
is deposited in the "central region" of the collision. Most of the baryons are carried by the 
receding nuclei (the fragmentation region) leaving this central region almost baryon free. 
Estimates of the energy density in this region give e ~ 3 — 5Gev/fm 3 corresponding to 
temperatures T « 200Mev. Immediately after the collision, hard scattering of quarks and 
gluons dominate the dynamics the gluons have mean-free paths estimated to be of order 
0.5fm and the quarks 1 — 2fm (the difference is mainly due to color factors) hence after a 
time of the order of about lfm/c the plasma is thermalized. 



The next stage of the dynamics is described by Bjorken's hydrodynamic picture |20 
When the plasma has achieved local thermodynamic equilibrium and for wavelengths longer 
than the mean free paths, the plasma can be described as a strongly coupled fluid and a 
hydrodynamic description is suitable. The essential ingredients in a hydrodynamics descrip- 
tion are i) the fluid is described by a local four velocity vector w M = r y(l,v);u^u^ l = 1, the 
energy momentum tensor is that of an homogeneous and isotropic fluid 
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= (e + P)u»u u - Pg^ 

with e , P the energy density and pressure respectively. The dynamics is then completely 
determined by conservation laws: i) baryon number, ii) energy momentum and by local 
thermodyamic equilibrium relations. The resulting picture of this hydrodynamic evolution 
is that the plasma expands and cools adiabatically and the temperature drops in time with 
the following law 

T(t)=T (^ 



T 

with c s the adiabatic sound speed, T > 200Mev and r ~ lfm/c. For a radiation dominated 
fluid c\ = 1/3. 

As the critical temperature for the chiral phase transition is reached from above the 
long-wavelength fluctuations of the chiral order parameter are expected to become critically 
slowed down. If this is the case the chiral phase transition can occur in a "quenched manner" 
and strongly out of equilibrium. Under these circumstances, Wilczek and Rajagopal argued 
that large domains in which the chiral order parameter could be "disoriented" with respect to 



the vacuum could appear |2T| • These domains are coherent pion condensates that form after 
the non-equilibrium phase transition much in the same manner as the correlated domains 
in condensed matter systems. These pion condensates decay, the neutral pion decays into 
two photons and the charged pions decay into muons. The pions can then be reconstructed 
and therefore these disoriented chiral condensates could lead to experimentally observable 
anomalies in the ratio of the number of neutral to charged pions. In isospin symmetric 
condensates the probability for finding a ratio R of neutral to total (neutral plus charged) 
is P(R) oc S(R — 1/3) (for large number of pions) whereas a disoriented chiral condensate 
leads to P(R) oc 1/y/R 0. 

The possibility of formation of disoriented chiral condensates had been previously con- 
jectured by Bjorken p|J as a potential explanation of CENTAURO events these are 



cosmic rays events with anomalous neutral to charged pion ratios. This possibility of a 
distinct signature associated with the chiral phase transition sparked an intense theoreti- 



cal effort [2^,p5[|. Several experimental searches are trying to find evidence for this pion 
condensates or disoriented chiral condensates at CERN-SPS (experiment WA98 ||26|| ) at 
the Tevatron at Fermilab (Minimax experiment [0]), the PHENIX and STAR detectors at 



RHIC [p8] at BNL can provide an event-by-event analysis of this potential candidates and 
the ALICE experiment scheduled at CERN-LHC includes the detector CASTOR ^ that 
will be studying CENTAURO type events. 

These disoriented chiral condensates are coherent pion domains and describe the same 
type of phenomenon as domains in quenched ferromagnets or in He superfluids as described 
previously. If these condensates are realized during a non-equilibrium stage of the chiral 
phase transition, they could lead to important experimental probes of this transition and 
hopefully will be amenable of detection at the RHIC and LHC on an event by event basis. 

There is an important difference in the dynamics of the chiral phase transition in the 
Early Universe and at Ultrarelativistic Heavy Ion Colliders. In the Early Universe, the 
chiral phase transition occured at a temperature of 150Mev when the Universe was about 
10 -5 seconds old in the radiation dominated era. The size of the Universe at that time was 
about 10 Km which is much larger than the mean-free path of quarks and gluons ~ 10 _15 m 
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and the time scale for cooling T/T pa 10~ 5 secs is much longer than the relaxation time 
scale of partons T re i pa 10~ 23 secs. Therefore in the Early Universe the confining and chiral 
phase transition most likely occured in equilibrium. These time and length scales must be 
compared to those in heavy ion collisions: the time scale for cooling from hydrodynamic 
expansion is few fm/c and the relaxation time scale near phase transitions could be longer 
and comparable to the lifetime of the quark-gluon plasma, furthermore current numerical 
estimates determine that the region in which the QGP is formed is about 20 fm. Hence there 
is a possibility that these phase transitions could be out of equilibrium in heavy ion collisions 
and that novel phenomena associated with the process of phase ordering and the emergence 
of pion condensates could be important experimental signatures of the chiral transition. 



C. Early Universe Cosmology: 

The COBE satellite mission revolutionized the field of Cosmology. The discovery of 
temperature fluctuations in the Cosmic Microwave Background (CMB) of 30fiK imprinted 
on a blackbody spectrum at 2.73K provides supporting evidence for the main ideas that 
seek to explain the small inhomogeneities that gave rise to large scale structure formation 
1 30jj3Tf1 . Two leading contenders with a solid base on particle physics have emerged that 



provide different explanations for the origin of the primordial inhomogeneities that grew via 
gravitational instability to form large scale structures: inflation and defects pOi^Ij fl. Infla- 



tion postulates that at an energy scale determined by Grand Unified Theories ~ 10 16 Gev 
the Universe underwent a period of exponentially accelerated expansion during which its 
size grew by a factor e 60 necessary to solve several problems with the standard Big Bang 
Cosmology p0|j31| . Small quantum fluctuations that were present during this epoch of infla- 



tion soon became causally disconnected and therefore unaffected by microphysical processes. 
These fluctuations became in causal contact again at a much later stage of the cosmological 
evolution, when the Universe was basically dominated by matter. Small fluctuations begin 
to grow under gravitational instability when they become causally connected again but after 
the epoch of radiation-matter equality at a temperature of about lOeV and redshift z pa 10 4 
|30|] . The COBE experiments are sensitive to those fluctuations that have established causal 
contact again after the epoch of recombination about 300000 years after the Big Bang at a 
redshift of about z ~ 1100. Therefore observations of the CMB allow to obtain information 
on the spectrum of primordial quantum fluctuations that were present shortly after the Big 
Bang. 



Original scenarios of inflation relied on a supercooled phase transition |30[. Recent 
detailed studies of the dynamics of phase transitions in early universe cosmology [32| allow 
a reliable calculation of the dynamics including backreaction effects on the metric and a 
self-consistent evolution of classical gravity and quantum fields. This approach allows to 
extract the power spectrum of the primordial perturbations of the metric arising from the 
quantum fluctuations in the matter fields. These fluctuations are directly related to those of 
the temperature of the CMB at the scale of recombination and correspond to the Sachs- Wolff 
plateau in the power spectrum measured by COBE [^,^l|]. It is found |32| that the growth 



of correlated domains after a supercooled phase transition of second order (no metastability) 



favors a power spectrum with more power on long wavelengths [32| as a consequence of the 



process of phase ordering. This is a consequence of the unstabilities associated with the 
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early stages of the phase ordering dynamics. This "red" power spectrum is consistent with 
the results of COBE for the temperature anisotropies, provided that the couplings of the 
matter field are fine tuned . 

The alternative proposal for explaining the small primordial fluctuations that result in 
the small temperature inhomogeneities observed on the CMB invoke the formation and 
subsequent evolution of a network of defects after a cosmological phase transition out of 
equilibrium P0|,5I|. In these scenarios, the dynamics of phase ordering after a phase transi- 



tion results in a network of "cosmic strings" |l]-§]]. The phase transitions required in these 
scenarios occur at a GUT scale ~ 10 16 Gev and the important quantity is the string tension 
H which is the energy per unit length of these topological defects and which determines 
the amplitude of the fluctuations on the space-time metric. There is a very important dif- 
ference between the fluctuations in the inflationary and the topological defects scenarios. 
In inflation, the quantum fluctuations become causally disconnected p0| , pl| and therefore 
their evolution is very simple until they become causally connected again because these 
fluctuations are not influenced by microphysical processes during the period of acausal evo- 
lution. Contrary to this "acausal" dynamics, topological defects are always causal and are 
constantly influenced by microphysical processes. Their evolution must be followed dynam- 
ically from the time at which the network of defects is formed, at a GUT scale, all the way 
up to the time scale at which they result in the formation of large scale structure-several 
billion years later!. Obviously this is an enormous dynamical range, however, detailed com- 
puter simulations reveal that a scaling solution emerges (for details see ||) determined by a 
dynamical length scale. The results of numerical studies suggest that this dynamical length 
scale is completely determined by the size of the causal horizon at a given time (see the later 
section on Cosmology for details on causal horizons). 

The emergence of this length scale through the dynamical process of phase ordering is 
exactly the same that has been previously discussed within the context of condensed matter 
systems. Current ground based and balloon borne experiments along with large scale surveys 
and future satellite missions will provide a flood of data that will support or falsify current 
theoretical ideas on large scale structure formation and temperature anisotropies. Thus an 
important theoretical effort goes in providing reliable predictions on the power spectrum 
of primordial quantum fluctuations. It is a tantalizing possibility that these cosmological 
observations could provide a definite evidence for cosmological phase transitions. 



D. ...Therefore 

We have seen in detail that the dynamics of phase ordering and evolution after non- 
equilibrium phase transitions are of fundamental importance in a wide range of energies from 
meV, in Condensed Matter, through Gev in the physics of the Quark Gluon Plasma and 
the Chiral Phase Transition all the way to GUT's (10 16 Gev) in Early Universe Cosmology. 
An important technical aspect in the study of these phenomena is their non-perturbative 
nature: in a rapid phase transition (of typical second order without metastable states) small 
amplitude long-wavelength fluctuations become unstable (this will be understood in detail 
below) and grow in time. The amplitude of these fluctuations must grow until they begin 
to sample the broken symmetry states of thermodynamic equilibrium. 
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Furthermore, in quantum field theory the notion of a "defect" requires a careful quan- 
titative understanding of the process of classicalization and decoherence of quantum fluc- 
tuations. Defects are intrinsically classical objects and are typically described as arising in 
a field theory as solutions of the classical equations of motion. Whether the description of 
the dynamics of quantum phase transitions in terms of the emergence of defects or coherent 
structures is a correct or appropriate one can only be determined by following the time evo- 
lution of the quantum density matrix. The classicalization of fluctuations and the emergence 
of semiclassical coherent, large amplitude field configurations should be a consequence of the 
time evolution and not an a priori description of the dynamics. Therefore one of our main 
goals is to provide a consistent quantum field theory description of the dynamics that is 
capable of describing this process of classicalization and decoherence of fluctuations. 

Having discussed in some detail the importance of the dynamics of symmetry breaking 
phase transitions out of equilibrium within important settings and their experimental study, 
we now provide some of the technical aspects that help clarify the phenomena and their 
quantitative study. We begin by describing a phenomenological approach to phase ordering 
kinetics in condensed matter systems, highlighting the important ingredients and concepts. 
We then move on to furnish a quantitative approach to the study of the non-equilibrium 
dynamics in quantum field theory to compare some striking similarities to condensed matter 
and also to contrast some important and relevant differences. The main point for delving 
into some technical details is to emphasize many robust features of the dynamics of symmetry 
breaking and phase ordering, 

• The early stages of phase ordering are determined by linear (spinodal) instabilities. 
Long-wavelength fluctuations become unstable and grow. 

• The emergence of a dynamical length scale. This scale represents the average size 
of the ordered domains and grows in time, eventually at asymptotically long times 
becoming macroscopically large. 

• Associated with this dynamical length scale there is dynamical scaling, asymptotically 
this length scale determines the behavior of correlation functions. 

• Coarsening: the growth in time of this correlation length translates in that the peak 
of the power spectrum moves towards longer wavelength, resulting in a sharp "Bragg 
peak" at asymptotically long times. This Bragg peak reflects the onset of condensates 
corresponding to ordered regions of macroscopic size. 

As we will see in detail, the phenomenological description in condensed matter systems is 
very different from the microscopic description in quantum field theory. Nevertheless we find 
that despite these important differences the above features are fairly robust and common to 
all of the situations studied. Only an excursion into the technical details can reveal in full 
force these very important and remarkable features. 

II. PHENOMENOLOGY OF PHASE ORDERING DYNAMICS IN CONDENSED 

MATTER: 

The phenomenological description of phase ordering kinetics begins with a coarse grained 
local free energy functional of a (coarse grained) local order parameter M(r) |7],§] which 
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determines the equilibrium states. In Ising-like systems this M(r) is the local magnetization 
(averaged over many lattice sites), in binary fluids or alloys it is the local concentration 
difference, in superconductors is the local gap, in superfluids is the condensate fraction etc. 
The typical free energy is (phenomenologically) of the Landau- Ginzburg form: 



F[M}= J d d x |^[VM(x)] 2 + V[M{x)\ 



V[M}= l -r(T)M 2 + ^M A ; r(T) = r (T - T c ) (2.1) 

Fig. 1 depicts V[M] for T > T c and T < T c . The equilibrium states for T < T c 
correspond to the broken symmetry states with M = ±Mq(T) with 

_ / _ for T >T C 

M {T)-\ /^ (Tc _ T) i forT<Tc (2-2) 



Below the critical temperature the potential V[M] features a non-convex region with 
d 2 V[M]/dM 2 < for 

- M S (T) < M < M 8 (T) ; M S (T) = ^(T - T c )^ (T < T c ) (2.3) 

this region is called the spinodal region and corresponds to thermodynamically unstable 
states. The lines M S (T) vs. T and M (T) vs. T [see eq.( j2.2|) 1 are known as the classi- 
cal spinodal and coexistence lines respectively. Fig. 2 displays the classical spinodal and 



coexistence curves for the potential V[M] in ( 2.1 



The states between the spinodal and coexistence lines are metastable (in mean-field 
theory). As the system is cooled below T c into the unstable region inside the spinodal, the 
equilibrium state of the system is a coexistence of phases separated by domains and the 
concentration of phases is determined by the Maxwell construction and the lever rule. 

Question: How to describe the dynamics of the phase transition and the process of 
phase separation? 

Answer: A phenomenological but experimentally succesful description, Time Dependent 
Ginzburg-Landau theory (TDGL) where the basic ingredient is Langevin dynamics 0- [Hj 

fl"M = _ r[ f,M] «3«L + , m (2.4) 
dt [ ' J 5M(r,t) /V ' ; V ; 

with r)(r, t) a stochastic noise term, which is typically assumed to be white (uncorrelated) 
and Gaussian and obeying the fluctuation-dissipation theorem: 

(v(r,t)ri(r',t')) = 2TT(r)5 3 (f-r )5(t - t') ; (r](f,t)) = (2.5) 

the averages (■ ■ •) are over the Gaussian distribution function of the noise. There are two 
important cases to distinguish: NCOP: Non-conserved order parameter, with T = To a 
constant independent of space, time and order parameter, and which can be absorbed in a 
rescaling of time. COP: Conserved order parameter with 

T[r\ = -T Vl 
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where r could depend on the order parameter, but here we will restrict the discussion to 
the case where it is a constant. In this latter case the average over the noise of the Langevin 
equation can be written conservation law 



dM 



Ot 



-V -J + T) 



d_ 

dt 



( / d 3 fM(r,t)) = 



J = V f 



-r 



5F[M) 
5M 



(2.6) 



where /i is recognized as the chemical potential. Examples of the NCOP are the magnetiza- 
tion in ferromagnets, the gap in superconductors and the condensate density in superfluids 
(the total particle number is conserved but not the condensate fraction), of the COP: the 
concentration difference in binary fluids or alloys. For a quench from T > T c deep into the 
low temperature phase T — > the thermal fluctuations are suppressed after the quench and 
the noise term is irrelevant. In this situation of experimental relevance of a deep quench the 
dynamics is now described by a deterministic equation of motion, 
for NCOP: 



dM 

~dT 



for COP: 



dM 

~~dT 



SF[M] 
5M 



5F[M] 
1 5M 



(2.7) 



(2.8) 



which is known as the Cahn-Hilliard equation []7|||. In both cases the equations of motion 
are purely diffusive 



dF 
~dl 



d 3 r 



SF[M] dM{r,t) 



8M(f, t) 



dt 



$d 3 r(gf NCOP 

/<Mvf£) 2 COP 



(2.9) 



and in both cases ^ < 0. Thus, the energy is always diminishing and there is no possibility 
of increasing the free energy. Thus overbarrier thermal activation cannot be described in 
the absence of thermal noise, which is clear since thermal activation is mediated by large 
thermal fluctuations. The fact that this phenomenological description is purely dissipative 
with an ever diminishing free energy is one of the fundamental differences with the quantum 
field theory description studied in the next sections. 



A. Critical slowing down in NCOP: 



Critical slowing down of long- wavelength fluctuations is built in the TDGL description. 
Consider the case of NCOP and linearize the TDGL equation above the critical temperature 
for small amplitude fluctuations near M = 0. Neglecting the noise term for the moment and 
taking the Fourier transform of the small amplitude fluctuations we find 



dm k (t) 
dt 



k 2 + r (T-T c ) m k (t) 



(2.10) 



12 



showing that long-wavelength small amplitud fluctuations relax to equilibrium m k 
time scale given by 



on a 



k 2 + r (T-T c ) 



(2.11) 



As T — > T 6 + the long-wavelength modes are critically slowed down and relax to equilibrium 
on very long time scales. Therefore a TDGL description leads to the conclusion that if the 
cooling rate is finite, the long-wavelength modes will fall out of LTE and become quenched. 
As the temperature falls below the critical, these modes will become unstable and will grow 
exponentially. 



B. Linear instability analysis: 



Let us consider now the situation for T «T C and neglect the thermal noise. The early 
time evolution after the quench is obtained by linearizing the TDGL equation around a 
homogeneous mean field solution M (t). Writing 



M(r, t) = M (t) + 4^E m k(t) e ifc> 



(2.12) 



where Q is the volume of the system, and considering only the linear term in the fluctuations 
rrik(t) the linearized dynamics is the following: COP: for M Q (t) the conservation gives 



dM (t) 
dt 







since M Q is the volume integral of the order parameter [see eq.( |2.6| )1 and for the fluctuations 
we obtain 



dm k {t) 
dt 



ui{k) mkit) ; oj{k) = —T k 2 



k 2 + 



d 2 V[M] 



dM 2 



M a 



(2.13) 



In the spinodal region 9 



d 2 V[M] 




dM 2 


M„ 


iield grow 



M 



< there is a band of unstable wave vectors k 2 < 



for which the frequencies are positive and the fluctuations away from the mean 



NCOP: separate the k ^ from the k = in the linearized equation of motion: 
dM Q (t) ^ dV[M) 



dt 

dm k (t) 
dt 



-r 
-r„ 



dM 

5F[M] 
5M 



m k {t) 



M (t) 



k 2 + 



d 2 V[M] 



dM 2 



Mo 



(2.14) 
(2.15) 



whereas the first equation ( |2.14| ) determines that M {t) rolls down the potential hill towards 
the equilibrium solution, the second equation also displays the linear instabilities for the same 
band of wave vectors as in the COP in the spinodal region |M (t)| < M S (T) [see eq. (|2.3| )1 for 
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which the fluctuations grow exponentially in time. Thus in the linearized approximation for 
both NCOP and the COP the spinodal instabilities are manifest as exponentially growing 
fluctuations. These instabilities are the hallmark of the process of phase separation and are 
the early time indications of the formation and growth of correlated regions which will be 
understood in an exactly solvable example below. 



C. The scaling hypothesis: dynamical length scales for ordering 

The process of ordering is described by the system developing ordered regions or domains 
that are separated by walls or other type of defects. The experimental probe to study 
the domain structure and the emergence of long range correlations is the equal time pair 
correlation function 

C(f,t) = (M(f,t)M(0,t)) (2.16) 

where (• • ■) stands for the statistical ensemble average in the initial state (or average over the 
noise in the initial state before the quench) and will become clear(er) below. It is convenient 
to expand the order parameter in Fourier components 



Ak-x 

k 

and to consider the spatial Fourier transform of the pair correlation function 

S(k,t) = (m^m^t)) (2.17) 

known as the structure factor or power spectrum which is experimentally measured by 
neutron (in ferromagnets) or light scattering (in binary fluids) Jl2|]. The scaling hypothesis 



introduces a dynamical length scale L(t) that describes the typical scale of a correlated 
region and proposes that 

t) = f {w) ) ^ t] = Ld{t) 9{kL{t)) (2 - 18) 

where d is the spatial dimensionality and / and g are scaling functions. Ultimately scaling 
is confirmed by experiments and numerical simulations and theoretically it emerges from a 
renormalization group approach to dynamical critical phenomena which provides a calcula- 
tional framework to extract the scaling functions and the deviations from scaling behavior [|7| . 
This scaling hypothesis describes the process of phase ordering as the formation of ordered 
'domains' or correlated regions of typical spatial size L(t). For NCOP typical growth laws 
arc L(t) « t 1 ' 2 (with some systems showing weak logarithmic corrections) and L(t) ~ t 1 ' 3 
for scalar and ~ t 1//4 for vector order parameter in the COP case 



D. An exactly solvable (and relevant) example: the Large N limit 

We consider the case where the order parameter has iV-components and transforms 
as a vector under rotations in an N-dimensional Euclidean space, i.e. M(f, t) = 
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(Mi(r, t), M2(r, t), ■ ■ ■ , Mjv(r, i)). For iV = 1 an example is the Ising model, for N = 2 
superfluids or superconductors (where the components are the real and imaginary part of 
the condensate fraction or the complex gap respectively), N = 3 is the spin one Heisenberg 
antiferromagnet, etc. For N = 1 the topological defects are domain walls (topological in 
one spatial dimension), for N = 2 they are vortices in d = 2 and vortex lines in d = 3, for 
N = d = 3 the topological defects are monopoles or skyrmions which are possible excita- 
tions in Quantum Hall systems and also appear in nematic liquid crystals 0. For iV — > oo 
and fixed d no topological defects exist. However the exact solution of the large N model 
gives insight and is in fairly good agreement with growth laws for fixed N systems which 
had been studied experimentally and numerically |?|,|T0|. In cosmological space-times it has 



been implemented to study the collapse of texture- like configurations H J33|J3^1 (see later). 
In quantum field theory the non-equilibrium dynamics of phase transitions has been studied 
in Minkowsky and cosmological space-times |25| , |35| -p9|] . The large N limit is an exactly 
solvable limit that serves as a testing ground for establishing the fundamental concepts and 
that can be systematically improved in a consistent 1/N expansion. It provides a consistent 
formulation which is non-perturbative, renormalizable and numerically implementable and 
has recently been invoked in novel studies of non-equilibrium dynamics in quantum spin 
glasses and disordered systems E0| . 



The exact solution for the dynamics in the large N limit, being available both in the 
condensed matter TDGL description of phase ordering kinetics and in Quantum Field The- 
ory in Minkowsky and Cosmological space times, allow us to compare directly the physics 
of phase ordering in these situations. Thus we begin by implementing this scheme in the 
NCOP case for the TDGL description. 

What is the (• • •) in the equations of the previous section?: consider that before the 
quench the system in in equilibrium in the disordered phase at T » T c and with a very 
short correlation length (£(T) ~ 1/T). The ensemble average in this initial state is therefore 

(M\r, 0)M j {7, 0)) = A 5 ij 5 3 {f - 7) 

(M i (f,0)) = (2.19) 

where A specifies the initial correlation. Now consider a critical quench where the system is 
rapidly cooled through the phase transition to almost zero temperature but in the absence 
of explicit symmetry breaking fields (for example a magnetic field). The average of the 
order parameter will remain zero through the process of spinodal decomposition and phase 
ordering. During the initial stages, linear instabilities will grow exponentially with m l k (i) ~ 
m l k (0) e^ fc )' ; u>(k) = k 2 — r(0) for k 2 < r(0) and at early times 

(mi(t)m^(t)) « A e 2 ^' (2.20) 

hence fluctuations begin to grow exponentially and eventually will sample the broken sym- 
metry states and the exponential growth must shut-off. The large N limit is implemented 
by writing the potential term in the free energy as 

r(T) - A -> -> 

V[M] = — ^M 2 + — (M 2 ) 2 ; M 2 = M-M (2.21) 

where A is kept finite in the large N limit. We will focus on the NCOP case with a quench 
to zero temperature and rescale the order parameter, time and space as 
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'r(O) 



M= V A ^ ; r (°) r ° t = r ' \M°) f=2; ( 2 - 22 ) 
after which the evolution equation for the NCOP case becomes 

where derivatives are now with respect to the rescaled variables. The large iV limit is solved 
by implementing a Hartree-like factorization f?j] 

r/ 2 — ► (r/ 2 ) = N(rjl) no sum over i (2.24) 

Then for each component the NCOP equation becomes 

V 2 + M 2 {t)]r H (2.25) 



drjj 
dr 

M\t) = 1 - ( V 2 ) (2.26) 



the eq. (|2.26|) is a self- consistent condition that must be solved simultaneously with the 
equation of motion for the components. Thus the large N approximation linearizes the 
problem at the expense of a self-consistent condition. The solution for each component is 
obviously 

vat?) = v a ( Q ) e ' k2T+b{T) ; Kr) = f m\t')<1t> (2.27) 



Consider for a moment that the k = mode is slightly displaced at the initial time, then it 
will roll down the potential hill to a final equilibrium position for which M 2 (oo) r)i(oo) = (so 
the time derivative vanishes in equilibrium). If 771(00) 7^ is a broken symmetry minimum of 
the free energy, then M 2 (r) — > when r — > 00. This is the statement of Goldstone's theorem 
that guarantees that the perpendicular fluctuations are soft modes. This asymptotic limit 
allows the solution of the self-consistent condition 

M 2 (r) = 1 - ( V 2 (r)) = 1 - A e 2b ^ [ e' k ' T = 1 - A e 2b ^ (87rtH (2.28) 

J (27r) d 

The vanishing of the right hand side in the asymptotic time regime leads to the self-consistent 
solution 



6(r)- j In 



M 2 (r) A (2.29) 



where ro is a constant related to A. This self-consistent solution results in the following 
asymptotic behavior 



d 



^)^^(°)(-J 6 ( 2 - 3 °) 

Introducing the dynamical length scale L{r) = t5 it is straightforward to find the structure 
factor and the pair correlation function 
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S(k,t) ocL d (t) e - 2(feLW)2 (2.31) 



C(f, t)(xe^) ; L(t)=t* (2.32) 



This behavior should not be interpreted as diffusion, because of the L d {t) in eqn. (|2.31 ) 
which is a result of the self-consistent condition. 
Important Features: 

• The 'effective squared mass' M 2 {t) 0: asymptotically there are massless excita- 
tions identified as Goldstone bosons. 

• Since M 2 (t) — > asymptotically, the self-consistent condition results in that (M 2 ) — > 
Nr(0)/yX, i.e. the fluctuations sample the broken symmetry states, which are equi- 
librium minima of the free energy. These fluctuations begin to grow exponentially at 
early times due to spinodal instabilities. 

• A dynamical correlation length emerges Lit) = t 1 ^ 2 which determines the size of the 
correlated regions or 'domains'. A scaling solution emerges asymptotically with the 
natural scale determined by the size of the ordered regions. These regions grow with 
this law until they become macroscopically large. Although this a result obtained in 
the large N limit, similar growth laws had been found for NCOP both analytically 
and numerically for N — 1 etc. M 



Coarsening: The expression for the structure factor (|2.31 ) shows that at large times 



only the very small wavevectors contribute to S(k,t), however the self-consistency 
condition forces the J ' k d ~ x dk S(k,t) — > constant thus asymptotically k d ~ x S(k,t) is 
peaked at wavevectors k « L~ x (t) with an amplitude L d (t) thus becoming a delta 
function S(k,t) 5 d (k). The position of the peak in S(k,t) moving towards longer 
wavelength is the phenomenon of coarsening and is observed via light scattering. At 
long times a zero momentum condensate is formed and a Bragg peak develops at 
zero momentum, this condensate however grows as a power of time and only becomes 
macroscopic at asymptotically large times. Coarsening is one of the experimental 
hallmarks of the process of phase ordering, revealed for example in light scattering 



1^] and is found numerically in many systems j7|. Thus the large N limit, although 
not being able to describe topological defects offers a very good description of the 
ordering dynamics. 



III. PHASE ORDERING IN QUANTUM FIELD THEORY I: MINKOWSKI 

SPACE-TIME 

A. A quench in Q.F.T. 

Although the phenomenological Time Dependent Landau Ginzburg theory is a succesful 
description of phase ordering kinetics in condensed matter systems, there is no first principle 
derivation from a microscopic theory of these equations of motion. Whereas microscopic 
descriptions either based on classical or quantum Hamiltonians lead to time reversal invariant 
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equations of motion, the TDGL equations are first order in the time derivative and therefore 
purely dissipative. 

A first principles, microscopic description of a quantum theory must begin with the 
Heisenberg equations of motion for operators or the Schroedinger or quantum Liouville 
equations for the quantum states or density matrix that describes the system. In this 
section we provide an introduction to the treatment of strongly out of equilibrium situations, 
in particular that of a "quench" in a quantum field theory system. 



This is the situation studied in JO]] for the dynamics of formation and evolution of 
disoriented chiral condensates during the chiral phase transition. 

The dynamics is completely determined by the microscopic field theoretical Hamiltonian. 
For a simple scalar theory the Hamiltonian operator is given by 

H = |rf 3 a;{in 2 (f,t) + i[V$(x,t)] 2 + y[$(f,t)]} (3.1) 

where $ is the quantum mechanical field and IT its canonical momentum. We want to 
describe a quenched scenario where the initial state of the system for t < is the ground 
state (or density matrix, see later) of a Hamiltonian for which the potential is convex for 
all values of the field, for example that of an harmonic oscillator, in which case the wave 
function(al) \E r [$] is a Gaussian centered at the origin. At t = the potential is changed so 
that for t > it allows for broken symmetry states. This can be achieved for example by 
the following form 

V[<S>] = l m 2 (t)$ 2 + ^$ 4 (3.2) 

2m / +m 2 > fort<0 , , 

m(t) = \-mi<0 fort>0 (3 ' 3) 

thus the potential in Fig. 1 changes suddenly from T > T c to T < T c . Although in 
Minkowski space-time this is an ad-hoc choice of a time dependent potential that mimics 



the quench ||42|| , we will see in the next section that in a cosmological setting the mass 
term naturally depends on time through the temperature dependence and that it changes 
sign below the critical temperature as the Universe cools off. Most of the results obtained 
in Minkowski space-time will translate onto analogous results in a Friedmann-Robertson- 
Walker cosmology. Unlike the phenomenological (but succesful) description of the dynamics 
in condensed matter systems, in a microscopic quantum theory the dynamics is completely 
determined by the Schrodinger equation for the time evolution of the wave function or 
alternatively the Liouville equation for the evolution of the density matrix in the case of 
mixed states. We will cast our study in terms of a density matrix in general, such a density 
matrix could describe pure or mixed states and obeys the quantum Liouville equation 

.dm 



i- 



H(t),p(t) (3.4) 



dt 

Question: How does the wave function(al) or the density matrix evolve after a quench? 
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B. A simple quantum mechanical picture: 



In order to gain insight into the above question, let us consider a simple case of one 
quantum mechanical degree of freedom q and the quench is described in terms of an harmonic 
oscillator with a time dependent frequency u: 2 {t) = —e(t) u 2 ; u 2 > with e(t) the sign 
function, so that u 2 (t < 0) > ; u 2 {t > 0) < 0. Furthermore let us focus on the evolution of 
a pure state (the density matrix is simple the product of the wave function and its complex 
conjugate). Consider that at t < the wave function corresponds to the ground state of the 
(upright) harmonic oscillator. For t > the wave function obeys 



M[q,t] 
1 dt 



Id 2 1 2 2 

W " 2"° q 



*[q,t] (3.5) 



Since the initial wave function is a gaussian and under time evolution with a quadratic 
Hamiltonian Gaussians remain Gaussians, the solution of this Schrodinger equation is given 
by 

= N(t) e~^ q2 (3.6) 

Separating the real and imaginary parts of A(t) it is straightforward to find that 
|iV(t)| 4 /Re[A(t)] is constant, a consequence of unitary time evolution. Eq.( |3.8| ) can be cast 
in a more familiar form by a simple substitution 

A(t) = =► 0(t) - uj 2 0(t) = (3.8) 



where the equation for was obtained by inserting the above expression for A(t) in fl3.8|) . 
The solution is 0(t) = a * + b e _a;o * featuring exponential growth. This is the quantum 
mechanical analog of the spinodal instabilities described in the previous section. The equal 
time two-point function is given by 

(q 2 )(t) = A R \t) = \m 2 ~e 2 ^ (3.9) 

The width of the Gaussian state increases in time (while the amplitude decreases to maintain 
a constant norm) and the quantum fluctuations grow exponentially. As the Gaussian wave 
function spreads out the probability for finding configurations with large amplitude of the 
coordinates increases. These is the quantum mechanical translation of the linear spinodal 
instabilities. When the non-linear contributions to the quantum mechanical potential are 
included the single particle quantum mechanical wave function will simply develop two peaks 
and eventually re-collapse by focusing near the origin undergoing oscillatory motion between 
'collapses' and 'revivals'. In the case of a full quantum field theory there are infinitely many 
degrees of freedom and the energy is transferred between many modes. This simple quantum 
mechanical example paves the way for understanding in a simple manner the main features 
of a quench in the large N limit in quantum field theory, to which we now turn our attention. 
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C. Back to the original question: Large N in Q.F.T. 



We now consider the large N limit of a full Q.F.T. in which 

<£(£, t) = , t), $ 2 (f , £),-••, $ N (x, t)) (3.10) 

and similarly for the canonical momenta II. The Hamiltonian operator is of the form ( ft.lp 
with 

V[$]= 1 -m 2 (t)$-$ + ^[$-$] 2 (3.11) 

with m 2 (t) given by (|3.3|) . Let us focus on the case in which the initial state pure and 
symmetric, i.e. (<&) = 0, with < • • • > being the expectation value in this initial state. The 
more complicated case of a mixed state, described by a density matrix is studied in detail 
in [j3"E|-|3~8| and the main features are the same as those revealed by the simpler scenario of a 
pure state. The large N limit is implemented in a similar manner as in the TDGL example, 
via a Hartree like factorization 

($■ $) 2 -> 2 ($■$)$■ $ (3.12) 

where the expectation value is in the time evolved quantum state (in the Schrodinger picture) 
or in the initial state of the Heisenberg operators (in the Heisenberg picture). Via this 
factorization the Hamiltonian becomes quadratic at the expense of a self-consistent condition 
as it will be seen below. It is convenient to introduce the spatial Fourier transform of the 
fields as 

$(x,t) = -j=£*if(*)e** (3.13) 

with Q the spatial volume, and a similar expansion for the canonical momentum U(x,t). 
The Hamiltonian becomes 

H = E {^n r n_ s + V fc 2 (t) * s ■ $_ fe -} (3.14) 

fe 

W 2 (t) = m 2 (t) + k 2 + — j ^ <$ r $_ 5 >(t) (3.15) 
The problem now has decoupled in a set of infinitely many harmonic oscillators, that are 



only coupled through the self-consistent condition in the frequencies (|3.15|) . To induce a 
quench, the time dependent mass term has the form proposed in eq. ( p.3|) . 

Just as in the simple quantum mechanical case, we consider the initial state to be a 
Gaussian centered at the origin in field space, which is the ground state of the (upright) 
harmonic oscillators for t < 0. Since a Gaussian is always a Gaussian under time evolution 
with a quadratic Hamiltonian, we propose the wave function(al) that describes the (pure) 
quantum mechanical state to be given by 

*[$,t] = u k hf k (t) <r^M-«} ; A k (t = o) = w k (t < o) (3.16) 
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Time evolution of this wavefunction(al) is determined by the Schrodinger equation: in the 
Schrodinger representation the canonical momentum becomes a differential (functional) op- 
erator, 11^ — > — i5/5§ and the Schrodinger equation becomes a functional differential 
equation. Comparing the powers of in this differential equation, one obtains the follow- 
ing evolution equations for N k (t) and A k (t) 

±]nN k (t) = - l -A k (t) (3.17) 

i^& = Al{t)-W2{t) (3.18) 

As in the single particle case, the constancy of \N k (t)\ 4 /Re[A k (t)} is a consequence of unitary 
time evolution. The non-linear equation for the kernel A k (t) can be simplified just as in the 
single particle case by writing 

<Pk{t) 

and taking the expectation value of $ 2 in this state we obtain 

($ s >$_ s )(t) = N\<l> k (t)\ 3 (3.20) 

Hence we find a self-consistent condition much like the one obtained in the large N limit for 
TDGL. The equations for the mode functions and the self-consistent condition for t > are 
therefore given by 

Mt) + [k 2 + M 2 (t)] <f> k (t) = (3.21) 
M 2 {t) = -ml + \j^\Mt)\ 2 (3-22) 

where the integral in the self-consistent term in (|3.22[) is simply ($1). There are two funda- 
mental differences between the quantum dynamics determined by the equations of motion 
and the classical dissipative dynamics of the TDGL phenomeno logical description given in 
sec. II: 

• The equations of motion and the self-consistency condition equations ( 3.21]) -( 3.22|) lead 
immediately to the conservation of energy 

• The evolution equations are time reversal invariant. 

These properties must be contrasted to the purely dissipative evolution dictated by the 
TDGL equations as is clear from eq. (2.9). Consider a very weakly coupled theory A << 



1 and very early times, then the self-consistent term can be neglected and we see that 
for k 2 < m\ the modes grow exponentially. This instability again is the manifestation of 
spinodal growth [ 43 - fl5| , j36| ,j37|| ■ Since the mode functions grow exponentially, fairly soon, at 



a time scale t s ~ mg 1 ln(l/A) the self-consistent term begins to cancel the negative mass 
squared and M 2 (t) becomes smaller. We find numerically that this effective mass vanishes 
asymptotically, as shown in Fig. 3. 
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D. Emergence of condensates and classicality: 



The physical mechanism here is similar to that in the classical TDGL, but in terms of 
quantum fluctuations. The quantum fluctuations with wave vectors inside the spinodally 
unstable band grow exponentially, these make the (<3> 2 ) self-consistent field to grow non- 
perturbatively large until when ($ 2 ) « tuq/X wnen the self-consistent (mean) field begins to 
be of the same order as m 2 , (the tree level mass term). At this point the quantum fluctuations 
become non-perturbatively large and sample field configurations near the equilibrium min- 
ima of the potential. The spinodal instabilities are shutting off since the effective squared 
mass M 2 (t) is vanishing. 

When M 2 (t) vanishes, the equations for the mode functions become those of a free 
massless field, with solutions of the form 4>k{t) = A^ e lkt + Bk e~ tkt , whereas for the k = 
mode the solution must be of the form <po(t) = a + bt with a; b ^ since the Wronskian of 
the mode function and its complex conjugate is a constant. This in turn determines that 
the low k (long wavelength) behavior of the mode functions is given by 

sin kt 

4>k\t) — a cos kt + b — - — (3.23) 

K 

This behavior at long wavelength has a remarkable consequence: at very long time the 
power spectrum \<fik{t)\ 2 , which is the equivalent of S(k,t) for TDGL (see eq. ( |2.17| )) is 



dominated by the small fc-region, in particular k « 1/t, with an amplitude that grows 
quadratically with time. Then the structure factor S(k,i) = \<fik(t)\ 2 features a peak that 
moves towards longer wavelengths at longer times and whose amplitude grows with time in 
such a way that asymptotically J °° k 2 S(k,t)dk/2n 2 — > m^/X and the integral is dominated 
by a very small region in k that gets narrower at longer times. This is the equivalent of 
coarsening in the TDGL solution in the large N limit, where the asymptotic time regime 
was dominated by the formation of a long-wavelength condensate. Fig. 4 shows the power 
spectrum at two (large) times displaying clearly the phenomenon of coarsening and the 
formation of a non-perturbative condensate. 



The pair correlation function can now be calculated using this power spectrum ]37 



C(f,t) = — — / ksmkr \4> 2 k (t)\ dk . (3.24) 
2 ix 1 r Jo 

At long times and distances the integral is dominated by the very long wavelength modes, 
in particular by the term oc sin[A;t]/A; of <f>k{t), hence the integral can be done analytically 
and we find 

C(r,t) = -G(2t-r) (3.25) 
r 

with A a constant. This is a remarkable result: the correlation falls off as 1/r inside domains 
that grow at the speed of light. This correlation function is shown in Fig. 5 at several 
different (large) times. This correlation function is of the scaling form: introducing the 



dynamical length scale L{t) = t it is clear that [37[ 



C{r,t)<xL-\t)f{r/L{t)) ; f( 8 ) = ^—^ (3.26) 
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We interpret these 'domains' as being a non-perturbative condensate of Goldstone 
bosons, with a non-perturbatively large number of them oc 1/A, such that the mean square 
root fluctuation of the field samples the (non-perturbative) equilibrium minima of the poten- 
tial. In particular an important conclusion of this analysis is that the long-wavelength modes 
acquire very large amplitudes, their phases vary slowly as a function of time (for k << 1/t), 
therefore these fluctuations which began their evolution as being quantum mechanical, now 
have become classical. 



E. O.K. ..O.K. but where are the defects? 

At this point our analysis begs this question. To understand the answer it is convenient 
to back track the analysis to the beginning. The initial quantum state is given by a the 
wave-function(al) (|3.16|) , thus the most probable field configurations found in this ensemble 
are those whose spatial Fourier transform are given by 

|$ fc | cc 1 ex 1 (3.27) 
'W k (t < 0) Jk 2 + ml 



(restoring h would multiply by y%). Then typical long- wavelength field configurations 
that are represented in the quantum ensemble described by this initial wave-function(al) are 
of rather small amplitude. The initial correlations are also rather short ranged on scales 
uiq . Under time evolution the probability distribution is given by 

( Kwi a ) 

P[$,t] = |^[<M| 2 = n^IIfc I \N k (t)\ 2 e W 1 \ (3.28) 

At times longer than the regime dominated by the exponential growth of the spinodally 
unstable modes, the power spectrum |0f.(t)| 2 obtains the largest support for long wavelengths 
k « 771q and with amplitudes ~ ttiq/X. Therefore field configurations with typical spatial 
Fourier transform 4>k{t) are very likely to be found in the ensemble. These field configurations 
are primarily made of long- wavelength modes and their amplitudes are non-perturbatively 
large, of the order of the amplitude of the fields in the broken symmetry minima. A typical 
such configuration can be written as 

&(x,t) typical « \Mt)\ co # ■ ^+ 41 ( 3 - 29 ) 

k 

where the phases 6% are randomly distributed with a Gaussian probability distribution since 
the density matrix is gaussian in this approximation. We note that a particular choice of 
these phases leads to a realization of a likely configuration in the ensemble that breaks trans- 
lational invariance. In fact translations can be absorbed by a change in the phases, thus 
averaging over these random phases restores translational invariance. Since the quantum 
state (or density matrix) is translational invariant a particular spatial profile for a field con- 
figuration corresponds to a particular representative of the ensemble. Combining all of the 
above results together we can present the following consistent interpretation of the ordering 
process and the formation of coherent non-perturbative structures during the dynamics of 
symmetry breaking in the large N limit 
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• The early time evolution occurs via the exponential growth of spinodally unstable 
long wavelength modes. This unstable growth leads to a rapid growth of fluctuations 
($ 2 )(t) which in turn increases the self-consistent contribution and tends to cancel the 
negative mass squared. The effective mass of the excitations -Wg + ^($ 2 )(t) — > 
and the asymptotic excitations are Goldstone bosons. 

• At times larger than the spinodal time t s w mo 1 ln(l/A), the effective mass vanishes 
and the power spectrum or structure factor S(k,t) = \4>k{t)\ 2 displays the features of 
coarsening: a peak that moves towards longer wavelengths and increases in amplitude, 
resulting in a long-wavelength condensate at asymptotically long times. 

• For large time a dynamical correlation length emerges L(t) = t and at long distances 
the pair correlation function is of the scaling form C(r,t) oc L~ l (t)f(r/L(t)). The 
length scale L(t) determines the size of the correlated regions and determines that 
these regions grow at the speed of light. Inside these regions there is a non-perturbative 
condensate of Goldstone bosons with a typical amplitude of the order of the value of 
the homogeneous field at the equilibrium broken symmetry minima. 

The similarity between these results and those of the more phenomenological TDGL descrip- 
tion in condensed matter systems is rather striking. The features that are determined by 



the structure of the quantum field theory are |37| : i) the scaling variable s = r/t with equal 



powers of distance and time is a consequence of the Lorentz invariance of the underlying 
theory, ii) the fact that the pair correlation function vanishes for r > 2t is manifestly a con- 
sequence of causality. An analysis of the correlations and defect density during the spinodal 
time scale has been performed in [H6l and related recent studies had been performed in p7| . 



IV. PHASE ORDERING IN QUANTUM FIELD THEORY II: FRW COSMOLOGY 

A. Cosmology 101 (the basics): 

On large scales > 100 Mpc the Universe appears to be homogeneous and isotropic as 
revealed by the isotropy and homogeneity of the cosmic microwave background and some of 
the recent large scale surveys Q. The cosmological principle leads to a simple form of the 
metric of space time, the Friedmann-Robertson- Walker (FRW) metric in terms of a scale 
factor that determines the Hubble flow and the curvature of spatial sections. Observations 
seem to favor a flat Universe for which the space time metric is rather simple: 

ds 2 = dt 2 -a 2 {t)dx 2 (4.1) 

the time and spatial variables t, x in the above metric are called comoving time and spatial 
distance respectively and have the interpretation of being the time and distance measured 
by an observer locally at rest with respect to the Hubble flow. At this point we must note 
that physical distances are given by l p h ys {t) — a {t) x. An important concept is that of 
causal (particle) horizons: events that cannot be connected by a light signal are causally 
disconnected. Since light travels on null geodesies ds 2 = the maximum physical distance 
that can be reached by a light signal at time t is given by 
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ft dt' 

d H {t) = a(t) / — - (4.2) 
Jo a(t') 

It will prove convenient to change coordinates to conformal time by defining a conformal 
time variable 

ft dt' 

rj= -f-^ds 2 = C 2 {r 1 ) {drf-dx 2 ) ■ C(rj) = a{t{rj)) (4.3) 



a(t') 



in terms of which the causal horizon is simply given by dsiji) = C{v) V an d physical 
distances as x p h ys = C(rj) x. This metric is of the same form as that of Minkowski space 
time. For energies well below the Planck scale Mpi ph 10 19 Gev gravitation is well described 
by classical General Relativity and the Einstein equations: 

BT - - g^R = (4.4) 

where we have been cavalier and set c = 1 (as well as h = 1). R^ v is the Ricci tensor, R 
the Ricci scalar and T^ v the matter field energy momentum tensor. The above equation 
is classical but one seeks to understand the dynamics of the Early Universe in terms of a 
quantum field theory that describes particle physics, thus the question: what is exactly the 
energy momentum tensor?, in Einstein's equations it is a classical object, but in QFT it is 
an operator. The answer to this question is: gravity is classical, fields are quantum mechan- 
ical, but T^ LU — > (T^), i.e. it is the expectation value of a quantum mechanical operator 
in a quantum mechanical state. This quantum mechanical state, either pure or mixed is 
described by a wave-function(al) or a density matrix whose time evolution is dictated by the 
quantum equations of motion: the Schrodinger equation for the wave functions or the quan- 
tum Liouville equation for a density matrix. Consistency with the postulate of homogeneity 
and isotropy requires that the expectation value of the energy momentum tensor must have 
the fluid form and in the rest frame of the fluid takes the form (T^ v ) = diagonal(p,p, p,p) 
with p the energy density and p the pressure. The time and spatial components of Einstein's 
equations lead to the Friedman equation 



d 2 (t) 8tt 



pi 



■Pit) (4-5) 



a 2 (t) 3M- 
d(t) a? it) 8tt 

2 W) + W) = -W p[,] 



PI 



Combining these two equations one arrives at a simple and intuitive equation which is 
reminiscent of the first law of thermodynamics: 

^(pa 3 (t)) = -p^2 p + 3% + p) = (4.7) 
dt dt a 

The alternative form shown on the right hand side of ( |4.7|) is the covariant conservation 
of energy. Since the physical volume of space is Vo a 3 (t) (with Vq the co moving volume) 
the above equation is recognized as dU = —p dV which is the first law of thermodynamics 
for adiabatic processes. To close the set of equations and obtain the dynamics we need 
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an equation of state p = p(p): two very relevant cases are: i) radiation dominated (RD) 
with p = p/3 and matter dominated (MD) p = (dust) Universes. In our study we will 
focus on the RD case. The equation of state for RD is that for blackbody radiation for 
which the entropy is S — CVT 3 (with C a constant). Since V(t) = V a 3 (t) is the physical 
volume, the equation (|4.7|) which dictates adiabatic (isoentropic) expansion leads to a time 
dependence of the temperature: T(t) = To/a(t). Now the cooling is done by the expansion 
of the Universe and a phase transition will occur when the Universe cools below the critical 
temperature for a given theory. For the GUT transition T c w 10 16 Gev « 10 29 i^, for the 
EW transition T c « 100 Gev « 10 15 K. Returning now back to the large N study of the 
dynamics of phase transitions, we can include the effect of cooling by the expansion of the 
Universe by replacing the time dependent mass term m 2 (t) in (|3.11|) by 



m 2 (t) = ttIq 



T 2 (t) 

c 



- 1 



This form is consistent with the Landau- Ginzburg description including the time depen- 
dence of the temperature via the isentropic expansion of the Universe, but perhaps more 
importantly it can be proven in a detailed manner from the self-consistent renormalization 
of the mass in an expanding Universe [|35| . Thus the large N limit in a RD FRW cosmology 
will be studied by using the potential ( p.llj ) but with the time dependent mass given by 



B. Large N in Radiation (RD) and Matter (MD) dominated FRW Cosmology 

The large N limit is again implemented via the Hartree-like factorization ( |3.12j ) perform- 
ing the spatial Fourier transforms of the fields and their canonical momenta and including 
the proper scale factors, the Hamiltonian now becomes [3l| 

H (t) = E {^y n s ■ n_ s + wlit) s s ■ (4.9) 



k 2 \ 



Wm = - + m\t) + — <* s ■ (4.10) 

where now the expectation value is in terms of a density matrix p[$Q, 3>(T); t] since we are 
considering the case of a thermal ensemble as the initial state. 

We propose the following Gaussian ansatz for the functional density matrix elements in 



the Schrodinger representation |35 



P [$, $, t] = n w) {-^ ^ • + ^ i s • i_ s + s s • *-^} (4- n ) 

This form of the density matrix is dictated by the hermiticity condition pt[$, $, £] = 
p*[$,$,t]; as a result of this, B^(t) is real. The kernel Bk(t) determines the amount of 
mixing in the density matrix, since if Bk = 0, the density matrix corresponds to a pure state 
because it is a wave functional times its complex conjugate. The kernels Afc(0) ; Bk(0) are 
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chosen such that the initial density matrix is thermal with a temperature T, > T c |35[ . Fol- 
lowing the same steps as in Minkowski space time, the time evolution of this density matrix 
can be found in terms of a set of mode functions <p k {t) that obey the following equations of 
motion and self-consistency condition 



(j> k {t) + 3 - fc (t) + 

(X 



k' 2 



a 2 {t) 



+ m 2 (t) 



m 2 (t) = ml 



T? 



+ 



A 



d 3 k 







(2k) 



\Mw coth 



2T, 



(4.12) 
(4.13) 



This equations can be cast in a more familiar form by changing coordinates to conformal 
time (see eq. (|4.3|) ) and (conformally) rescaling the mode functions <f>k{t) = fk{xi)/C(rj) 
obtaining the following equations for the conformal time mode functions fk(v) i n a FRW 
cosmology 



f'M+\k 2 + C 2 ( V )M 2 (r ] )}f k ( V )=0 



M 2 (eta) = ml 



T 2 



1 



d 3 k 



\Mv)f coth w k (o) 



c 2 (v) 



2T, 



C"{rj) 



(4.14) 
(4.15) 



where primes now refer to derivatives with respect to conformal time. For RD and MD FRW 

V 



C(rj) 
C{rj) 



1 + 7T ! C"(r]) = for RD 



:i + 



C"(rj) = 1/8 for MD 



(4.16) 

(4.17) 
(4.18) 



where we have rescaled the only dimensionful variable m , length and time are now in terms 
of wig 1 . The above equations of motion now have a form analogous to those in the case of 
Minkowski space-time. 

As the temperature falls below the critical the effective squared mass term becomes 
negative and spinodal instabilities trigger the process of phase ordering. This results in that 
the quantum fluctuations quantified by ($ 2 ) grow exponentially. These spinodal instabilities 
make the self-consistent field grow at early times and tends to overcome the negative sign 
of the squared mass, eventually reaching an asymptotic regime in which the total effective 
mass M 2 (rf) vanishes. 

Again this behavior determines that the fluctuations are sampling the equilibrium broken 
symmetry minima of the initial potential, i.e. ($ 2 ) — > 2N ™° . 

Although, just as in Minkowski space-time the effective mass vanishes asymptotically, the 
non-equilibrium evolution is rather different. We find numerically |39| that asymptotically 
the effective mass term behaves as 



C 2 { V )M 2 {ri) -Ih/Arj 2 for RD 
C 2 (t])M 2 (t]) -35/4?7 2 for MD 



(4.19) 
(4.20) 



Fig. 6 displays C 2 (f])M 2 (r]) as a function of conformal time for the case of Tj/T c = 1.1 
with T c oc mo/VX US for RD. 
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We see that at very early time the mass is positive, reflecting the fact that the initial 
state is in equilibrium at an initial temperature larger than the critical. As time evolves the 
temperature is red-shifted and cools and at some point the phase transition occurs, when 
the mass vanishes and becomes negative. 

Figure 7 displays 2I } 2 (& 2 )(v) vs - V i* 1 units of m^ 1 for |f = 3, g — 1CT 5 for an R.D. 

c 

Universe. Clearly at large times the non-equilibrium fluctuations probe the broken symmetry 
states. 

This particular asymptotic behavior of the mass determines that the mode functions 
fk(v) grow as rf/ 2 for RD and for rfl 2 k < \jr\ and oscillate in the form e ±tkri for k > 1/r). 
This behavior is confirmed numerically []3"5|1 . We find both analytically and numerically f3"5fl 
that asymptotically the mode functions are of the following form in terms of the scaling 
variable x = kr) 

f k ( V ) = A V I ^ + B^H 2 (x) for RD (4.21) 

f k ( V ) = Arfi ^ + B^H 3 (x) for MD (4.22) 

x A 772 

Where A,B are numerical constants and ^2,3(2;) are Bessel functions. For fixed x, and 
r] » 1 the asymptotic behavior is completely determined by the J 2 ,3(x), i.e., asymptotically 
the solutions are of the scaling form. 

Figure 8 displays f]~ 5 \fk{v)\ 2 as a function of the scaling variable krj revealing the scaling 



behavior in RD, a similar behavior emerges for MD [p9| . The absolute value of the coefficient 
A is completely fixed by the sum rule 

A (&)(V) = 1 (4-23) 



2Nm 2 

in particular this sum rule determines that \A\ 2 oc 1/g. 

It is remarkable that this is exactly the same scaling solution found in the classical non- 
linear sigma model in the large N limit and that describes the collapse of textures [33], and 
also within the context of TDGL equations in the large N limit applied to cosmology |3~3 . 
However, there are important differences that are purely quantum mechanical in origin: the 
Bessel functions H 2 ^(x) in the solutions (|4.21| , fii~22j ) are necessary to maintain the constancy 



of the Wronskian for the mode functions, the coefficients A, B contain information of the 
initial state. The contribution of the if 2,3(2?) breaks the scaling property of the solutions, but 
they are necessary for a consistent solution for the effective mass term that falls off as v^ 2 . 
In particular this behavior of the effective mass constrains the non-perturbative contribution 
to B of 0(l/y/g) [|39J. 

The growth of the long-wavelength modes and the oscillatory behavior of the short 
wavelength modes again results in that the peak of the structure factor S(k,r]) = \ fk{v)\ 2 (X 
C 2 (f])r] 3 g(kri) moves towards longer wavelengths and the maximum amplitude increases. 
This is the equivalent of coarsening and the onset of a condensate. 

Although quantitatively different from Minkowsky space time, the qualitative features 
are similar. Asymptotically the non-equilibrium dynamics results in the formation of a non- 
perturbative condensate of long-wavelength Goldstone bosons. We can now compute the 
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pair correlation function C(r, rj) from the mode functions solutions to ( |4.14 ) and find that it 



is cutoff by causality at r = 2r\. The correlation function computed with the mode functions 
in the asymptotic regime agrees perfectly with that computed from the asymptotic form 
given by (|4.21 ). The correlation function is depicted in Fig. 9 for two different (conformal) 



times. 

The scaling form of the pair correlation function is 

C(r,rj) oc rf x( r /2v) 

where xi x ) is a hump-shaped function as shown in fig. 9. 

Clearly a dynamical length scale L{rj) = rj emerges as a consequence of causality, much 
in the same manner as in Minkowsky space time. The physical dynamical correlation length 
is therefore given by Cphysiv) = C{rj)L{rj) = d H (t), that is the correlated domains grow again 
at the speed of light and their size is given by the causal horizon. The interpretation of this 
phenomenon is that within one causal horizon there is one correlated domain, inside which 
the mean square root fluctuation of the field is approximately the value of the equilibrium 
minima of the tree level potential, this is clearly consistent with Kibble's original observation 



IJfJ]. Inside this domain there is a non-perturbative condensate of Goldstone bosons |39 
There are remarkable consequences of the scaling solution 



When the scale factor C(rj) oc T] a (a = 1, 2 for RD and MD respectively) the asymptotic 
solution is universal in the sense that it does not depend on the past history of the 
background metric, the coefficient A is completely determined by the sum rule an is 
insensitive to the past evolution of the scale factor. This is important in the case of 
a transition between RD and MD, if the two regimes are well separated in conformal 
time, the asymptotic solution in the RD and MD eras is independent of the transition 



• We find |39[ that the fluid that results from the fluctuations obeys the same equation 
of state of the background fluid. In particular we find p = | — lj, i.e. p = e/3 for 
a = 1(RD) and p = for a = 2 (MD). 

Thus we have seen that the phenomenon of scaling, coarsening and the onset of conden- 
sates during the non-equilibrium dynamics of phase ordering is a universal feature of the 
process of phase ordering. The non-perturbative large N limit has allowed a clear comparison 
between the phenomenological description in condensed matter based on the TDGL, and 
the microscopic quantum field theoretical description in Minkowski and FRW space-times. 



V. CONCLUSIONS AND LOOKING AHEAD 



In this lectures we have discussed the multidisciplinary nature of the problem of phase 
ordering kinetics and non-equilibrium aspects of symmetry breaking. Main ideas from con- 
densed matter were discussed and presented in a simple but hopefully illuminating frame- 
work and applied to the rather different realm of phase transitions in quantum field theory 
as needed to understand cosmology and particle physics. In particular we have emphasized 
robust features of the process of phase ordering kinetics: early stages dominated by spinodal 



29 



instabilities and the growth of correlated regions, the emergence of a dynamical correlation 
length that determines the size of the correlated regions as a function of time and dynamical 
scaling at long times. The phenomenon of coarsening is a result of this scaling behavior and 
is reflected in that the peak in the power spectrum moves towards longer wavelengths, and 
asymptotically long times results in a "Bragg peak" that signals the onset of macroscopic 
ordered phases and condensates. The study of condensed matter systems was in terms of 
the phenomenologically succesful Time Dependent Landau Ginzburg theory which is purely 
dissipative and for which there is no first principles derivation from a microscopic theory in 
general. 

We then passed onto the study of the dynamical evolution out of equilibrium in quantum 
field theories both in Minkowsky and FRW space-times by providing a consistent non- 
perturbative framework to study the time evolution of an initially prepared density matrix. 

The large N approximation has provided a bridge that allows to cross from one field to 
another and borrow many of the ideas that had been tested both theoretically and exper- 
imentally in condensed matter physics. There are, however, major differences between the 
condensed matter and particle physics-cosmology applications that require a very careful 
treatment of the quantum field theory that cannot be replaced by simple arguments. The 
large N approximation in field theory provides a robust, consistent non-perturbative frame- 
work that allows the study of phase ordering kinetics and dynamics of symmetry breaking 
in a controlled and consistently implementable framework, it is renormalizable, respects all 
symmetries and can be improved in a well defined manner. This scheme extracts cleanly 
the non-perturbative behavior, the quantum to classical transition and allows to quantify 
in a well defined manner the emergence of classical stochastic behavior arising from non- 
perturbative physics. The emergence of scaling and a dynamical correlation length, coarsen- 
ing and the onset of non-perturbative condensates are robust features of the dynamics and 
the Kibble-Zurek scenario describes fairly well the general features of the dynamics, albeit 
the details require careful study, both analytically and numerically. 

We have emphasized that this study has very definite potential experimental implications, 
in QCD if the chiral phase transition occurs out of equilibrium in ultrarelativistic heavy ion 
collisions leads to the possibility of formation of disoriented chiral condensates that are 
described in the same manner as ordering domains in condensed matter. These condensates 
have a very distinct hallmark in that they lead to a very different ratio of neutral to charged 
pions, this property can be measured on an event by event basis with the detectors at the 
forthcoming heavy ion colliders. 

In cosmology the process of formation of ordered regions that grow after a rapid phase 
transition, the emergence of scaling and a dynamical length scale and coarsening of these 
domains lead to a definite prediction of a "red" power spectrum on scales that have re- 
enterd the causal horizon right after recombination. These are the scales that contribute to 
the temperature anisotropies measured by COBE and the forthcoming cosmological exper- 
iments. Therefore the study of the dynamics of symmetry breaking out of equilibrium in 
quantum field theory directly bears on experimental possibilities in a wide range of energies 
both in accelerator and cosmological experiments and is therefore an endeavour that must 
be pursued vigorously. 

Of course this is just the beginning, we expect a wealth of important phenomena to be 
revealed beyond the large N, such as the approach to equilibrium, the emergence of other 
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time scales associated with a hydrodynamic description of the evolution at late times and 
a more careful understanding of the reheating process and its influence on cosmological 
observables. Although within very few years the wealth of observational data will provide 
a more clear picture of the cosmological fluctuations, it is clear that the program that 
pursues a fundamental understanding of the underlying physical mechanisms will continue 
seeking to provide a consistent microscopic description of the dynamics of particle physics 
and cosmological phase transitions. 
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FIG. 5. grC{r,r) vs r/\m \ for t/\m \ = 100, 200, 300 for g = 10~ 7 . 
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FIG. 7. 2jVTO 2 (^H 7 ?) vs - 77 (conformal time in units of ttiq 1 ) for ^ = 3, g = 10 5 . R.D. 
Universe. 
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FIG. 9. C(r,7]) vs. r for rj = 250,400 (in units of m^ 1 ) for J = 3, g = 10~ 5 . R.D. FEW 
Universe. 
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